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PREFACE 


The ODIN procedure is a design analysis technique which allows 
the use of existing computer codes as part of a larger simula- 
tion. Communication of information among computer codes is 
accomplished by means of a data base repository accessible and 
managed by the ODIN executive computer code, ODINEX. 

The objective of the contract was the development of improved 
independent geometry display program. The result is a system 
of computer codes for editing geometry monitoring geometric 
perturbations, and drawing realistic pictures of complex 
vehicle geometries. 
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NOTE ON MEASUREMENT SYSTEMS 


Several vehicle drawings contained in this report include dimensional 
and performance data to illustrate features of the Aerophysics 
Research Corporation graphics codes. These numbers are quoted in 
U. S. Customary Units. These numbers may be converted to metric 
equivalent units using the following factors. 


Quantity 

U. S. Unit 

Metric Unit 

Conversion 

Length 

feet 

meter 

0.30480 

Area 

square feet 

square meter 

0.092903 

Volume 

cubic foot 

cubic meter 

0.028317 

Weight 

pound 

kilogram 

0.453 

Thrust force 

pound 

kilogram 

0.453 

Specific 

impulse 

second 

second 

1.0 

Velocity 

feet/second 

meters/second 

0.30480 


It should also be noted that certain of the computer graphics program 
input quantities which serve to position a drawing on the plot paper 
are set to nominal values in inch units. These values must be quoted 
in inch units when using the codes. For those interested in converting 
metric units to U. S. Customary Units, the multiplicative conversion 
factor for centimeters to inches is 0.39370. 



SUMMARY 


The PCSYS programs use a vehicle geometric definition based 
upon quadrilateral surface elements to produce realistic 
pictures of an aerospace vehicle. It is anticipated that the 
PCSYS system will be an important component of the Optimal 
Design Integration (ODIN) system, see sketch below. The PCSYS 
programs can be used to visually check geometric data input, 
monitor geometric perturbations, and to visualize the complex 
spatial inter-relationships between the internal and external 
vehicle components. The pictures constructed with PCSYS can be 
annotated with text information if desired. 

PCSYS has two major component programs. The first program, 

IMAGE, draws a complex aerospace vehicle pictorial representation 
based on either an approximate but rapid hidden line algorithm 
or without any hidden line algorithm. The second program, 

HIDDEN, draws a vehicle representation using an accurate but 
time consuming hidden line algorithm. Program IMAGE is based 
on the widely distributed code of Gentry, ref. 1. HIDDEN is 
based on a code developed by researchers at Purdue University. 



SKETCH OF THE ODIN SYSTEM FOR VEHICLE DESIGN SYNTHESIS 


The contents of this report are arranged in three parts. Part I 
describes the original ODIN picture drawing program IMAGE. The 
second part describes the new picture drawing program HIDDEN as 
it has been used in the ODIN system. These adaptations consist 
of 

a. An increase in the number of panels which may be 
employed in the definition of a vehicle surface 

b. Conversion of the IBM 360 version of the program 
to the CDC 6600/7600 series computer 

c. Extension of the code to configurations of the 
complexity presented in this report. Primarily 
this extension consisted of debugging the existing 
code. 

Part III covers adaptation of the HIDDEN code to the ODIN system. 
Two programs have been written to accomplish this adaptation. 
These programs are CNVTHD and PLOTHD. Program CNVTHD 
automatically transforms the geometry input of the existing 
IMAGE code into the form accepted by HIDDEN. Program PLOTHD 
converts the output of HIDDEN into a form acceptable to the 
ODIN picture drawing programs. As a result of the present study 
effort, program HIDDEN has been completely integrated into the 
ODIN system. 
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PART I 


PROGRAM IMAGE 


The generation of geometry for aerospace configurations in 
digital format is one of the more tedious tasks in the 
design analysis process and is also difficult to check. The 
program described in Part I of this report provides a pictorial 
representation of the digital geometric input. This program 
capability provides a visual check for errors in geometry and 
a picture presentation which can be annotated with text infor- 
mation. 

The IMAGE program is extracted from the program of reference 1. 

It has undergone extensive modification and simplification in 
an effort to provide a more useful addition to the ODIN (Optimal 
Design Integration) program library, references 2 and 3. Al- 
though developed primarily for use in the ODIN system, the 
IMAGE program is equally useful as an independent program. The 
program was written for the CDC 6600 computer. Some coding 
peculiar to the CDC machine must be changed for use on other 
machines but this coding is a small portion of the total program. 

Versions of the IMAGE program are currently available for the 
CDC 6600 computer, reference 4, and the UNIVAC 1108, reference 5. 
The code is also operational on the CDC 7600 as a result of the 
present contract. Mr. Alan Wilhite of NASA's Langley Research 
Center has also prepared a version which operates on the "Prime" 
mini-computer . 
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SURFACE MODEL 


The surface shape of the configuration is described to the 
IMAGE program by ordered sets of points in three dimensional 
space. The geometry is specified in 80 column (BCD) card 
format. It may be input directly or generated by another 
program and passed to the IMAGE program as a binary or BCD 
file. The program rotates the geometric data to prespecified 
viewing angles then transforms it into a plane coincident with 
the plane of the paper. The geometry is reordered into quadri- 
lateral elements for plotting. 

A grouping of four surface points is used to describe a quadri- 
lateral surface element. An organization of a large number of 
related elements forms a component. A number of components may 
be used to give a complete description of the configuration. 

Each component is an independent unit of geometry which may be 
drawn separately or collectively with other components . In 
general, the equations for the geometry rotations described 
here apply to any orientation angle. The equations required 
to produce the perspective drawings are derived in the follow- 
ing paragraphs. Figure 1 is a collection of component geometries 
representing a shuttle orbiter configuration arranged in a three 
view drawing by the program. 

Coordinate System 

Each point on the surface is described by its coordinates in 
the body reference coordinate system. 


X 

Y 

Z 

The body reference coordinate system is assumed to be a con^ 
ventional right-handed Cartesian system as illustrated below: 



ODIN/IMflGE 

ORBITER DESIGN ?SK P/L 


PICTURE DRAWING PROGRAM 



VEHICLE CHRRfiCTERISTICS 


1. HRSS PROPERTIES 


LRNOED HEIGHT 
ENTRY HEIGHT 

design cg. subsonic. CSK P/L 

DESIGN CG HYPERSONIC. 13K P/L 


I5S93S .5 
I9BG1C.7 
69.S597J 

GS-HSCn 


2. GEOMETRY 


BODY LENGTH 
TOTAL HING AREA 
CHORD. THEO ROOT 
CHORD. TIP 
ASPECT RATIO 
LEADING EDGE SHEEP ANGLE 
'RAII ING EDGE SHl.I.'P ANaC 
El.EVON AREA. TOTAL 
FHPOSEa NING LOCATION 


IIO-OCCO 
0500. CCO 
SO.SSoSS 
6.8Sii‘l 
3 .o:g::3 
RS.c;::o 
0.0CCG03 
375.0030 
58.45033 


3. AEROOYNAmp CHARACTER I ST I C.S 


DESIGN TRIH LIFT COEFFICIENT .6257300 
DESIGN nIHIHUH LANDING SPEED I87.8S37 
HRH TRIH ALPHA HYPERSONIC DESIGN. CONO 88.63321' 



FIGURE 1 ILLUSTRATION OF PICTURE AND TEXT OPTIONS IN IMAGE 


Coordinate Transformations 

To create the perspective drawings illustrated in this report 
each surface point on the body must be rotated to the desired 
viewing angle and then transformed into a coordinate system in 
the plane of the paper. With zero rotation angles the body 
coordinate system is coincident with the fixed system in the 
plane of the paper. 


Z 

o 



The rotations of the body and its coordinate system to give a 
desired viewing angle are specified by a yaw-pitch-roll sequence 
(’l',0,(t)). The rotation is given by the following relationship: 



This sequence is important to remember when describing the 
desired viewing angles to the IMAGE program. 


6 


The rotation matrices rp, 6 and <{) are given by; 


or 


whore 




cosi/r sini/f o 
~ sini/f cosi/f o 
o o 1 


W = 

cos 6 o 
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sin 6 o 

- sin0 
o 

COS0 
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1 o 

O COS(|l 

o -sin<J> 

o 

sin<j> 

COS(|> 

X 



rx„i 

Y 

Z 

■ • 


M 

“■ 1 
o o 
N 

j 

H 

= 

H H [+ 



Since each point on the surface is given by its coordinates in 
the X, Y, Z system, its position in the fixed coordinate system 
(Xq, Yq, Z^) may be found by the inverse of the above process: 
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If this operation is carried out, the resulting relationship is 
obtained. 


rx 1 

o 


cos6cosvl< -sinijj cos<j>+sin9cos4' sin(j> sirnli sin4>+sin6cos4' cos4> 


-X' 

Y 

o 

= 

cos0simi» cosi^ cos<j>+8in6sinijj sin4> -cosi]; sin4>+sin0sin4/ cos4) 


Y 

,Z 

o 


-sin0 cos0sin4> cos0cosi^ 


.z. 


X -X(cosGcoi\i/) + Y(-3in4jcos4>+sin6cosUi3ir.c-} + Z (sirxusino-f sinGcos'l'COsc) 
o 

-X(co‘39sinV) t Y(cosv|(Coso+sir.O sin'^JSlno) ■l-Z{-cos^i/£.i•n•trrsinGsir.vJ/cos6) 
o 

=^X(-sin0) Y(cosGsin6) + Z(cosGcodo) 


We may now use these last two equations to transform a given 
point on the body (X, Y, Z) with a specified set of rotation 
angles (\p,6,(l>) into the plane of the paper (the Y, Z system). 

With the CALCOMP library subroutines it is a simple matter to 
plot these data and to connect the related points with straight 
lines . 

The above relationships completely describe the transforma- 
tion required to rotate every point on the vehicle to the desired 
angle, then into the plane of the paper. However, the resulting 
drawing is difficult to interpret because all hidden lines are 
drawn. Further, since only one half, or one quarter (for 
symmetrical vehicles) of the coordinate points are usually 
input, some additional calculations are desirable. 

The surface points are grouped into quadrilateral surface elements 
for analysis by the program. Each quadrilateral is drawn as 
an individual 'curve.' The collection of all 'curves' represents 
the complete configuration. The technique also provides a con- 
venient means of limiting the drawn lines to those normally seen 
by the viewer at the defined viewing angles. 

Each input element is replaced by a plane quadrilateral surface 
made up of the four lines connecting the points. The quadri- 
lateral characteristics are used to determine the visibility 
of the four lines. The quadrilateral characteristics include 
the area, centroid and the direction cosines of the surface 
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unit normal. The surface unit normals may be transformed 
through the required rotation angles just as was done for the 
individual points . The resulting value of the component of 
the unit normal in the X direction (out of the plane of the 
paper) may be found from® the following equation: 


nxo “ (cos e cos '|/)+n^(-sinijJcos<b+ sine cosipsin4>)+n^{sin4jsin<j3+sinecos4jcos4>) 


where n^, n^, n^ are the components of the surface unit normal 
in the vehicle reference system. 


If n^ is positive, then the surface element is facing the 
o 

viewer. If n is negative, the element faces away from the 
^o 

plane of the paper. This result is used in the program to pro- 
vide the option of deleting most of those elements on a vehicle 
that normally could not be seen by a viewer. The picture is 
thus made more realistic and easier to interpret. Confusing 
elements which are on the back side of a component do not 
appear.. No criterion is provided, however, for the deletion 
of those elements that face the viewer but are blocked by 
other body components. If desired, more realistic drawings 
may be obtained by selective deletion of sections or half 
sections and by a proper selection of viewing angle. It should 
be noted that in certain open bodies portions of the outline of 
a body may be defined by quadrilateral elements which face away 
from the viewer. In such cases portions of the body outline 
will be deleted by the approximate hidden line algorithm contained 
contained in IMAGE. This problem can be avoided by use of the 
true hidden line algorithm described in Part II of the present 
report. 
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PROGRAM USAGE 

The computer program usage requirements described in this 
section are generally oriented toward the GDC 6600 computer 
version and specifically toward the Langley Research Center 
(LRC) installation. The actual program input requirements 
described are largely applicable wherever the program is 
installed but the control cards for the retrieval and execution 
of the program will differ from computer to computer and from 
installation to installation. 

The use of the program requires three types of input; the con- 
trol cards, the actual program input and the post processing 
instructions for generating CALCOMP plots. Figure 2 illustrates 
the deck setup for executing the program and indicates the 
separate submittal of a plot request card. One input case is 
illustrated. Multiple cases may be computed by repeating 
the case data illustrated. 

Control Cards 

The program input data is preceded by the operating system con- 
trol cards required to retrieve from permanent storage and 
execute the program. Figure 3 illustrates the control cards 
for three different ways of using IMAGE at LRC. Figure 3A shows 
the method of execution from the stored machine language pro- 
gram. Figure 3B illustrates a compile, load and execute 
sequence from stored source code, assuming program modifications 
are desired. Figure 3C illustrates the use of IMAGE within the 
ODIN (Optimal Design Integration) system. See reference 2 for 
ODIN system usage. 

The dashes on the JOB, USER and REQUEST cards indicate missing 
information which the user must supply in accordance with LRC 
computer complex accounting procedures . Further , the wedge 
number (program storage location in data cell) is subject to 
change. The latest wedge number is available from LRC. 

The user must supply the estimated run time in CPU seconds, 
the octal field length for the job and the number of operating 
systems (0/S) calls. These are given in the above order on 
the JOB card. Typical values for a single case are tabulated 
below; 
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PLOT TERMINATOR 



FIGURE 2 


DECK SETUP FOR IMAGE. 




RUNID,T20, CM2 0000. 

ACCOUNT CARD 

CHARGE CARD 

GET , IMAGE/UN= 7 2 7 8 5 ON . 

IMAGE. 

PLOT. 

7-8-9 

(IMAGE CASE DATA) 
6-7-8-9 


FIGURE 3A EXECUTION OF STORED PROGRAM 


RUNID,T20,CM45000. 

ACCOUNT CARD 

CHARGE CARD 

FETCH , A3 6 4 7 , SOURCE . 

RUN,S, , ,SCFILE. 

LGO. 

PLOT. 

7-8-9 

(MODS TO SOURCE PROGRAM, IF ANY) 
7-8-9 

(IMAGE CASE DATA) 

6-7-8-9 


FIGURE 3B COMPILE, LOAD, AND EXECUTE 


'EXECUTE IMAGE' 
(IMAGE CASE DATA) 
7-8-9 

'EXECUTE PLOTSV 
7-8-9 


REQUIRED ONLY IF PLOTS 
ARE REQUESTED 


FIGURE 3C EXECUTE IMAGE WITHIN 


FIGURE 3. ILLUSTRATIONS OF CONTROL 

PROGRAM IMAGE 


ODIN SIMULATION 


CARDS REQUIRED FOR 
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CPU 


FIELD LENGTH 


EXECUTE ABSOLUTE 
BINARY PROGRAM 

20 

35000 

COMPILE, LOAD 
AND EXECUTE 

40 

45000 

EXECUTE WITHIN 
THE ODIN SYSTEM 

30 

56000* 


*Minimum size for the executive system (reference 2) . 

The tabulated values above will vary with the complexity of 
the configuration (number of elements) , the number of cases 
and the number of plots requested. A good rule to follow is 
to allow ample CPU and 0/S calls in the first run, then use 
the day file results to estimate these parameters for similar 
configurations . 


General Program Input 

The IMAGE program input will typically consist of four types 
of input. 

1. Title card in free field format. 

2. Program controls ($TYPE32) in NAMELIST format. 

3. Surface model data as formatted corner-point geometry. 

4. Picture specifications ($TYP3435) in NAMELIST format. 

The program controls and picture specifications are input in 
standard NAMELIST input format. Two NAMELIST reads are provided 
for this purpose as illustrated in figure 4. 

The use of NAMELIST input was chosen for the following reasons : 

1. It is a simple name oriented input easily understood 
by most computer users. 

2. The format is standard and does not require relearn- 
ing from program to program. 

3. It is easily modified by the engineer or programmer 
when adding input variables to the program. 
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I TITLE (one card - first 59 characters) 


$TYPE32 

j^Geometry Scaling and Control Options J 


’Element Data - only if unit 5 
was specified by the data set 
. above 


1 


$TYP3435 

j^Picture Drav^ing or Text Options 

$ 


$TYP3435 

[ Picture Drawing or Text Options" 

(if 

LAST=1 


I END OF IMAGE DATA 


col. 71 ►‘9 9 


FIGURE 4 ILLUSTRATION OF INPUT STREAM TO 
IMAGE FOR TWO VIEWS. 
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When a NAMELIST read is encountered in the program, the entire 
input file is scanned up to an end-of-file or a record with a 
dollar ($) in column 2 followed immediately by the NAMELIST 
name requested by the programs. Succeeding data items are 
read until a second dollar ($) is encountered signifying the 
end of the NAMELIST. Any data on the input file before the 
requested NAMELIST is found will be ignored. All data between 
the opening and closing dollar is interpreted by the NAMELIST 
input routine. The data item within the NAMELIST statement 
may be in any of three forms: 

V = C, 

A = 

A(n) = , 

V is a variable name; C is a constant; A is an array name and 
n is an integer constant subscript, D^,...D^, are simple con- 
stants or repeated constants of ' the form k*C, where k is the 
repetition factor. Constants may be real, integer, hollerith 
or logical. Hollerith constants are preceded by nH where n 
is the number of characters in the hollerith constant. Logical 
constants are of the form .TRUE, (or T) or .FALSE, (or F) . 

Data items and constants must be separated by commas. The 
number of constants, including repetitions given for an 
unscripted array must equal the number of elements in that 
array. For a subscripted array name, the number of constants 
need not be equal but may not exceed the number of array 
elements needed to fill the array. More than one card may be 
used for input data and arrays may be split between cards. 

All except the last record must end with a constant followed 
by a comma and no sequence numbers may appear. The first 
column of each record is ignored. The set of data items may 
consist of any subset of the variable names associated with 
the NAMELIST name and the name need not be any particular 
order. More details on the use of NAMELIST are available in 
any FORTRAN user's guide, but the above description should be 
sufficient for the operation of the IMAGE program. 

The first list ($TYPE32) is used for specifying general data 
related to orientation, translation and scaling of the geometry. 
It also provides for setting flags pertaining to printing and 
geometry file control. 
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The $TYPE32 data set is followed by element geometry data if 
(described below) the INPUT file is specified in the above 
data set. Alternately the data may be read from the internal 
unit number 8. The geometry data is read, translated and 
scaled according to user instructions, then placed on a scratch 
unit (TAPE3) in binary format for use in the remainder of the 
program. 

After geometry data (if any) the second NAMELIST ($TYPE3435) 
is input. This data set provides input for picture control 
options, etc., required to generate the desired pictures. 

Figure 4 illustrates an input stream to IMAGE. Any number of 
views may be specified by repeating the $TYP3435 data set. 

The last data set should specify LAST = 1. 

The parameter, LAST, terminates the picture sequence for the 
current geometry which was temporarily placed on TAPE 3. The 
program logic returns to read a new TITLE card. The input 
flow is illustrated in figure 5. If a TITLE card is present 
in the input stream, additional geometry will be read from 
TAPES or TAPES and placed on the temporary storage file, TAPE 3 . 
A new sequence of pictures will be expected after the geometry 
is read. The execution is terminated by placement of a special 
(type 99) card in the input stream in place of a TITLE card. 

The special card must contain the integer, 99 in columns 71 
and 72. The following paragraphs describe each input type in 
detail . 

Title Card. - The TITLE card must be the first card in the 
input sequence. This card may contain from 1 to 59 characters 
(columns) of information which will be printed at the top of 
the frame, 9-inches above and 3-inches to the right of the 
initial reference point. A TITLE card must be present for 
each set of geometry which is to be processed by the program. 
Failure to supply this card will result in an input error. If 
no title is desired, a blank card must be inserted. The 
characters (99) in columns 71 and 72 of the TITLE card will 
cause normal termination of the program. 

Program Controls . - The $TYPE32 data set is a NAMELIST input 
set consisting of input instructions for the geometry data, 
scaling, translation and orientation of the data, and printing 
instructions for the quadrilateral element characteristics. 
Figure 6 summarized the NAMELIST names and descriptions for 
the input set. 

This section discusses each NAMELIST input in detail. Each 
paragraph is headed by the name and default value for the 
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START 


READ TITLE CARD 
59 CHARACTERS 


% Print program header 
at the top, 9.5 inches 
above the reference 
point. 

% Title card is printed 
9 inches above the 
initial reference point 


TEST FOR 
"TYPE 99 



READ $ TYPE 3 2 
NAMELIST DATA 


* Options of scaling the 
vehicle data and trans- 
lation of the vehicle 
coordinate system. 


READ GEOMETRY 
FROM TAPES OR 
TAPES (ITAPE OPTION) 


% Geometry sections are 
scaled and merged into 
a single section and 
placed on TAPE 3. 


READ $TYP3435 
NAMELIST DATA 


% Instructions for scaling 
and plotting the picture 
data. 


LAST 


• Test for more geometry. 


FIGURE 5 ILLUSTRATION OF INPUT FLOW 
LOGIC FOR IMAGE. 





NAMELIST 

NAME 

DEFAULT 

VALUE 

DESCRIPTION 

( 

DELX 

1. 0 

X- translation of scaled input data. 

i 

1 

f 

DELY 

1. 0 

Y-translation of scaled input data. 

S 

DELZ 

1.0 

Z-translation of scaled input data. 

1 

1 

lORIEN 

0 

Integer definition of a element 
orientation. 

:i 

1 

j 



= 0 Cross section input mode. 

1 



= 1 Streamwise input mode. 

( 

a; 



= 2 or 3 See text. 

S 

1 

IREW8 

0 

Control integer for logical unit 8. 
file position. 

1 



= 0 Rewind tape 8 before reading. 
=1 Do not rewind 8 before reading. 

1 

ISTAT3 

1 

Number of vehicle components. 

\i' 

'.i 

Ji: 

ITAPE 

0 

Geometry control integer. 



= 0 Geometry from unit 5. 

= 1 Geometry from unit 8 (coded) . 

■ \ 
:l 
! 



= 2 Geometry from unit 8 (binary) . 

1 

1,' 

:!' 

PRINTS 

0 

Print flag. 

i 



=0 No printing. 

= 1 Print quadrilateral data. 

H'J 

1! 

ill 

s 

'j 

it 

XSC 

0. 0 

X-scale factor. 

ii 

‘1 

YSC 

0. 0 

Y-scale factor. 

l! 

ZSC 

o 

• 

O 

Z-scale factor. 


FIGURE 6 

$TYPE32 NAFIELIST INPUT SET. 

1 
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variable listed. The heading is followed by a description of 
the input variable. If the default value is acceptable, the 
user need not define a value for it in the $TYPE32 NAMELIST 
data set. 

PRINTS = 0 PRINTS is an integer variable which may have mean- 
ingful values of zero (0) or one (1) . The variable 
controls the printing of detailed quadrilateral 
element characteristics. 

0 Element characteristics will not be printed. 

1 Element characteristics will be printed. 

lORIEN = 0 lORIEN is an integer variable which defines the 
element orientation for the data. All vehicle 
components must have the same orientation. 

0 Normal mode using cross sections. 

1 Geometry is input in streamwise strips. 

2 Geometry is input in streamwise strips. For each 

strip of elements, the first coordinate point in 
the right-hand strip of points is not used in the 
formation of the leading edge element but is 
ignored by the program. 

3 Same as = 2 except the left-hand point is ignored 
in the formation of the leading edge elements. 

Usually the data is described with lORIEN values of 0 or 1. 

Values of 2 or 3 will give correct pictures when the data-point 

slip methods of reference 1 are used to input the geometry data. 
Vehicle components with different values for lORIEN cannot be 
drawn correctly on a single picture frame. 

The scale factors and associated translation increments described 
below are generally used to transform the configuration geometry 
to a more convenient form for plotting. The factors are 
frequently used to move the vehicle reference axis so that the 
vehicle center corresponds approximately to the coordinate 
system origin. For a vehicle with its nose at X = 0.0 this 
is accomplished by using a DELX value (see below) of about one- 
half of the vehicle length. This simplifies the selection of 
picture scales to keep the vehicle within the imaginary picture 
frame . 
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The original geometry data on tape 5 or tape 8 is not changed 
by the use of the scale factors. They are applied to the 
geometry stored on the temporary file TAPES . The data is 
transfoinned using the following equations for X, Y and Z. 


X 

new 

Y 

new 

Z 

new 

The above 
XSC =1.0 

YSC = 1.0 

ZSC = 1.0 

DELX = 0.0 

DELY =0.0 

DELZ =0.0 


^input * 

(XSC) 

+ 

DELX 

^input ■ 

(YSC) 

+ 

DELY 

Z . ^ . 

input 

(ZSC) 

+ 

DELZ 


scaling terms are defined as follows: 

Scale factor to be multiplied by 

Scale factor to be multiplied by 

Scale factor to be multiplied by 

X increment to be added to X. . (XSC). 

input 

Y increment to be added to Y^j^^put ’ • 

Z increment to be added to Z . , . (ZSC). 

input 


Usually the configuration geometry consists of more than one 
component. In the data format of the IMAGE program, the data 
for each component is terminated with an integer flag referred 
to as a status flag as described under Surface Model Data 
(below) . The merging of the geometric components is desirable 
for plotting purposes. Therefore, the geometric components 
are merged into a single component as the geometric data is 
being transferred to the temporary file, TAPES. The user of 
the program must specify the number of components to be merged 
for each picture sequence. This specification is accomplished 
by the integer input variable ISTAT3. 

ISTAT3 = 1 Integer variable number of vehicle components 
with Status = 3 in the vehicle geometry which 
the user wishes to plot as a unit. The pro- 
gram will count the number of Status = 3 in 
the geometry deck and when the count reaches 
this input value, the program will proceed to 
the plot options for the vehicle components 
which have just been read. 
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Several options are available to the user of the program for 
accessing geometric data. Alternate files may be employed and 
alternate formats may be specified. The input parameter for 
controlling the above options is the integer variable ITAPE. 

ITAPE = 0 Geometry tape control integer variables with 
the following possible values: 

= 0 Geometry data (type 3) will be read from 
Tape 5 (geometry data cards are loaded 
along with picture-data control cards) . 

= 1 Geometry data (Type 3 will be read from the 
geometry storage tape (Tape 8) in coded 
format . 

=2 Geometry data (Type 3) will be read from 
the geometry storage tape (Tape 8) in 
binary format. 

The IMAGE program provides a flexible means of controlling the 
alternate geometry file, TAPE 8 . Multiple components may be 
stored on TAPES. These components may be extracted in sequen- 
tial groups or plotted individually. The input parameter for 
rewinding the geometry tape (TAPES) is IREW8 . Usually the 
file is rewound for the first sequence of pictures, then 
through the use of the input variable ISTAT3 , additional groups 
of components can be extracted for plotting purposes . 

IREW8 = 0 Integer variable to control the position of 

Tape 8 just before the geometry data are read 
from it. 

= 0 Rewind Tape 8 and then read geometry data 
from it. 

=1 Do not rewind Tape 8 , but start reading 

geometry data from it in its current position. 

The integer , IREW8 , permits the user to store more than one 
group of vehicle geometries on the geometry tape, then collect 
and plot them by groups . 

Surface Model Data . - The program accepts the geometry in 
element data form only. No simplified geometry techniques are 
employed. Several auxiliary programs are available which con- 
vert simplified geometry to element data format suitable for 
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input to the IMAGE program. Reference 1 is an example of such 
a program. It generates element data using several mathematical 
surface generation options and includes an aircraft geometry 
option which provides a convenient means of generating aero- 
dynamic surfaces. 

The element data in IMAGE may be read from input or from an 
alternate unit. If the INPUT file is selected, the actual 
data cards are merged with the other input data as illustrated 
in figures 2 and 4. 

If the alternate input unit (logical unit 8) is selected for 
element data input, the data may be read in coded (same as 
INPUT cards) or binary (fast read) mode. In order to use the 
binary mode, the data must have been written in binary mode. 
Reference 1 generates element data in coded format. The binary 
mode is very useful for improving efficiency particularily 
when the geometry must be regenerated many times for a problem 
solution. Regeneration of geometry is often desirable when 
studying the effects of geometric shapes. Regeneration will 
be required when optimizing geometric shape in an ODIN simula- 
tion. 

Element Data Format . - The element data method uses a large 
number of surface coordinate points on the surface of the 
configuration. The points can be ordered around a station 
contour or in a streamwise manner. Each point consists of 
an X, Y, Z coordinate set and a status flag. Each card con- 
tains two points. 

The coordinate system used for all the geometry data is shown 
in the figure below. For symmetrical vehicles it is standard 
practice to input the left side of the vehicle only. There 
is only one input method for the corner-point geometry and 
since any of the auxiliary geometry programs finally produce 
geometry data in surface-element form, it is important that 
the methods and nomenclature used with this method be clearly 
understood. It is, therefore, recommended that the input 
instructions for the surface-element method be studied before 
an attempt is made to use the method or write a geometry gen- 
eration routine. 
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The geometric input data in this method include the coordinates 
of a large number of points on the vehicle surface. The input 
data are organized in a manner that permits the description of 
a vehicle on a component buildup basis. This gives increased 
flexibility in shape description and makes it possible to draw 
exploded views by physically separating the components of the 
vehicle. Because of possible changes in the surface contours 
of a configuration, or replacement of the entire component, it 
may also be desirable to divide the configuration into several 
components. This permits easy changes either in a manual or 
automated mode such as ODIN. Each component of a configuration 
is further divided into a number of sections each defined by 
a group of points in space. In practice, the surface coordin- 
ates are usually recorded from cross-section drawings of the 
vehicle in such a way that each point need be read only once 
(even though it may be a member of as many as four adjacent 
quadrilateral elements) . Each point is defined by its three 
coordinates and a status flag that indicates whether it is 
the first point of a new section, a continuation of a group 
of points defining a section, the beginning of a new section, 
or the last point of the component. The program uses the 
status flags to determine how the input points are to be 
related to form the quadrilateral elements , and how the 
elements are combined to form a component. 

The first question that the user asks when starting to load the 
element geometry is the order in which the surface points are 
entered. The basic rules to be followed are given below. The 
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rules discussion will be followed by a discussion of a visual 
technique that many users will find helpful in determining the 
proper loading order. 

For the purpose of organizing the input data for computation, 
each point is assigned a pair of integers, m and n. These 
integers are not actually input to the program (they are cal- 
culated internally) but their use in the following discussion 
will provide a better understanding of the input data organiza- 
tion. For each point, n identifies the "column" of points to 
which it belongs, and m identifies its position in the "column," 
i.e., the "row." The first point of a "column" always has m = 1. 
To insure that the program will compute outward normal vectors, 
the following condition for the order to input points must be 
satisfied. If an observer is located in the outside the com- 
ponent and is oriented so that locally he sees points on the 
surface with m values increasing upward, he must also see n 
values increasing toward the right. Strict adherence to this 
simple rule will always lead to a correct set of input geometry 
data. Examples of correct and incorrect input are shown in 
the sketches below. In these pictures the exterior of the con- 
figuration lies above the paper, and the interior of the con- 
figuration lies below the paper. The arrows indicate the order 
of reading the points . 
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Associated with each input point is an input quantity called 
its status. The first point of each new section has Status 
= 2. Except for the first n-line of a component, the first 
point of each n-line has Status 1. The last point of the 
component of the vehicle has Status 3. All other. points have 
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Status = 0 (i.e., they may be left blank on the input sheet). 

The IMAGE program plots the picture according to components 
ending with a Status = 3 . 

The simple visual technique described below is helpful in 
determining the proper order of the input points : 

1. First, assume that you are holding in your hand a 
small model of the vehicle shape. Many program users 
find it helpful to construct a small paper model to 
help in visualizing the geometry loading procedure. 

On this model draw lines to represent the elements to 
be loaded for a given vehicle section. 

2. Next, decide which strips of elements are to constitute 
"columns" and which "rows." In most problems one of 
two procedures is selected - either a "column" of 
elements starts at the bottom of the shape and continues 
around to the top, roughly following vehicle cross- 
section lines, or a "column" is oriented so that it 
starts at the front part of the vehicle and runs aft 
toward the rear . 

3 . Hold the model out in front of you and rotate it until 
the columns are vertical with the first row of elements 
at the bottom. This procedure should be used regard- 
less of what part of the vehicle is being loaded - the 
body, fin, inside of fin, etc. Always orientate the 
model so that you are looking at the section to be 
loaded (from the outside, looking at the surface) with 
the columns running vertical , and the rows running 
horizontal . 

4 . Now that you have the section being loaded oriented 
in front of you, with the columns vertical, apply the 
following cardinal geometry rule: 

If a column of data points are loaded from the 
bottom to the top , then the next column of points 
(starting with a Status = 1) must be to the right. 

All of the geometry input data for this geometry option are 
input on "Type 3" element data cards (an integer 3 in column 
72). Each card contains the X, Y, Z coordinates and status flag 
for two points on the body surface. Every card in the element 
geometry deck must contain two surface points except the last 
card, which may have only the first surface point coordinates 
and status filled in. If a particular line of vehicle points 
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is odd in nximber, then it is usually advisable to repeat the 
last point (a dummy point) so that the last card will have 
two sets of point data. This permits the shifting of con- 
figuration components without disrupting other components. 

A description of the input card for element data is shown in 
figure 7 . 

Picture Specifications . - The picture control data set is a 
NAMELIST input set called $TYP3435 which is summarized in 
figure 8 . It controls the drawing of pictures and the print- 
ing of text. A typical deck setup will consist of several data 
sets of this type, one for each picture desired. The program 
will always try to read one data set. Additional views or text 
information may be generated or added by the inclusion of 
additional $TYP3435 data sets. The input integer controlling 
this function is LAST. 

LAST = 0 Integer variable controlling the program 

flow logic after drawing a picture or print- 
ing text. 

= 0 Return for $TYP3435 data set. 

= 1 Return for new TITLE card (may result 
in program termination or the reading 
of more geometric data) . 

If the variable LAST is set to one (1) , the reference axes 
system of the plotting device is moved to a new frame position 
specified by the user. The input variables controlling the axis 
translation are XMOVE and YMOVE. 

XMOVE = 17.0 The translation of the reference axis system 
in the X-direction (in inches) following 
the last $TYP3435 data set. 

YMOVE = 0.0 The translation of the reference axis system 
in the Y-direction (in inches) following the 
last $TYP343_5 data set. 

The values of XMOVE and YMOVE are made with respect to the 
original coordinate system reference specified upon execution 
of the program. The picture positioning parameters DXG and DYG 
described below have no affect on the XMOVE, YMOVE translation. 


The program permits the user to specify the viewing angles with 
the three-axis system discussed in the early part of this report 
The righthand rule is used for identifying the positive rotation 
angle of the vehicle geometry. The sequence of rotation is yaw, 
pitch and roll. 
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Column Code 


Explanation 


1-10 X 


X-coordinate of surface point (the 
value of X is written anywhere in 
this space with a decimal point and 
sign; usually input only if it is 
negative) . 


11-20 Y 


y-coordinate of surface point. 


21-30 Z 


Z-coordinate of surface point. 


31 STAT Status flag for the above set of 

coordinates (=2, 1, 0, or 3). 


32-41 


X-coordinate of surface point. 


42-51 


Y-coordinate of surface point. 


52-61 


Z-coordinate of surface point. 


STATT Status flag for the above set of 
coordinates (=2, 1, 0, or 3). 


66-68 CASE 


Case numl'jer (right- justified integer) 


69-70 SECT 


Numbers or letters to identify the 
vehicle section. These must be legal 
machine characters . 


72 TYPE 


Card type number = 3. 


FIGURE 7 ELEMENT DATA INPUT CARDS 
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NAMELIST 

NAME 

DXG 

DYG 

HTEXT 

lAREA 


ICS 


IQUAD 


I REEL 


ISHAD 


LAST 


FIGURE 8A 


DEFAULT 

VALUE DESCRIPTION 


0 . 0 
0.0 


Repositioning of X- and Y- reference 
point before plotting current 
picture “ inches. 


0.14 Character height for text material - in. 

0. Print control integer. 

= 0 No print. 

= 1 Print area of each section. 

0. Connectivity flag for quadrilaterals. 

= 0 Connect all four points. 

= 1 Connect points 1-2 and 3-4. 

= 2 Connect points 1-4 and 2-3. 

=4 Do not connect points. 

0. Controls points to be drawn. 

= 0 Draw input points. 

= 1 Draw computed points on 
quadrilateral . 

1. Reflected element flag. 

= 0 Draw input elements only. 

= 1 Draw reflected elements. 

= 2 Draw reflected elements (only 
one quadrant is input) . 

0. Hidden line option. 

=0 Do not plot hidden lines. 

= 1 Plot all lines. 

0. Program control flag. 

= 0 Return for $TYP3435 data set. 

= 1 Return for TITLE card (or 
Type 99 to terminate) . 


$TYP3A35 NAMELIST INPUT SET. 
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NAMELIST 

NAME 

PHI 

PSI 

THETA 


SCAL 

TEXT 


DEFAULT 

VALUE 

0 , 

0 . 

0 . 


DESCRIPTION 

Roll angle (see below) - degrees. 
Yaw angle (see below) - degrees. 
Pitch angle (see below) - degrees. 


Cl" yaw 


i,’,' 

}i_ 

l'!.'- 


// W: 


^0 > 




<t>, roll 


- 0 - 
0, pitch 


6.0 Frame size in inches. Geometry will 

be scaled to this dimension . 

F Logical control variable for text input. 

= .TRUE. Text information will follow 
$TYP3435. 

= .FALSE. No text information will be 
input. 


FIGURE 8B $TYP3435 NAMELIST INPUT SET. (CONTINUED) 
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NAMELIST DEFAULT 


NAME 

VALUE 

DESCRIPTION 

XLG 

* 

Extreme forward X-dimension on 
input geometry - input units. 

XRG 

¥ 

Extreme aft X-dimension on input 
geometry - input units. 

YBG 

* 

Extreme bottom Y-dimension on input 
geometry - input units. 

YTG 

* 

Extreme top Y-dimension on input 
geometry - input units. 

XMOVE 

17.0 

Translation of X-reference after 
one complete case - inches. 

YMOVE 

o 

. 

o 

Translation of Y-reference after 
one complete case - inches. 


* Computed 
read but 

before. This input set is 
may be overridden (see text) . 


FIGURE 8C $TYP3il35 NAMELIST INPUT SET (CONTINUED). 
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PSI = 0 . Yaw angle in degrees measured with respect 

to the positive Z-axis of the reference 
geometry coordinate system. 

THETA = 0. Pitch angle in degrees measured with respect 
to the positive Y-axis of the reference 
geometry coordinate system. 

PHI = 0. Roll angle in degrees measured with respect 

to the positive X-axis of the reference 
geometry coordinate system. Note that the 
positive X-axis is usually toward the nose 
of the vehicle. 

The technique used for plotting geometric data is to treat each 
quadrilateral element as a five-point 'curve' in the image 
plane. The IMAGE program provides the user with the option of 
specifying the points on each quadrilateral element which will 
be connected when plotting the elements. The integer variable 
which controls this option is ICS. 

ICS = 0 Integer variable controlling the points to 

be connected. 


= 0 Connect all 4-points of each element. 

= 1 Connect points 1-2 and 3-4 (see diagram 
below) 

= 2 Connect points 1-4 and 2-3. 

=3 Do not connect points with lines 



Quadrilateral 

Element 


Usually the geometric data plotted by the IMAGE program is 
symmetrical about the X-Z plane and data is generally provided 
for only one side of the configuration. The program provides 
the user the option of plotting only the input geometry, the 
reflected geometry or alternately plotting both the input 
geometry and the geometric reflection of the input geometry. 
The input integer controlling this option is IREFL. 

IREFL = 1 Integer variable controlling the plotting of 
reflection-elements as follows : 
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=0 Do not plot elements reflected to 
negative side of Y-axis. 

= 1 Plot reflected elements . 

= 2 Plot reflected elements (only one 
quadrant is input) . 

The plotting of all quadrilateral elements can often be distrac- 
ting if not misleading with respect to appearance of the 
vehicle which the geometric data represents. The distraction 
is usually caused by the plotting of elements which face away 
from the viewer. The IMAGE program provides the user with the 
option of eliminating the above class of elements for the 
individual geometry components. The integer variable which 
controls this option is ISHAD. IMAGE does not provide any 
capability for omitting elements which face the viewer but 
are hidden by another component of the vehicle. 

ISHAD = 0 Integer variable controlling the plotting of 

elements facing away from the viewer. 

=0 Do not plot elements that face away from 
the viewer (shadow elements) . 

= 1 Plot shadow elements (elements facing 
away from viewer) . 

The IMAGE program provides the user the option of printing the 
surface area characteristics of the vehicle components as follows: 

lAREA = 0 Integer variable controlling the printing of 

section areas. 

=0 Do not print section areas. 

= 1 Print out the area of each section. 

Usually, the desired picture is represented by the actual corner 
points of the geometric data. However, for some types of 
analysis, it is desirable to view the corner points of the 
actual quadrilateral elements since the latter represents the 
data actually used by some of the technology programs which the 
IMAGE geometry supports. The IQUAD option permits the selection 
of the points to be drawn. 

IQUAD = 0 Integer variable which selects the corner 
points to be drawn. 
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= 0 Draw input elements . 

= 1 Draw picture using quadrilateral element 
corner points. 

The program provides a flexible means of framing the picture 
data. Because of the flexibility provided, the user should 
exercise a certain degree of caution in setting up the data to 
assure the resulting picture will be plotted in the desired 
location. The following discussion of input variables is 
designed to help the unfamiliar user in setting up the data for 
positioning the picture sequences. Figure 9 shows the relation- 
ships among the variable input data discussed below. 

SCAL = 6. Imaginary frame size in inches for the current 
picture- All data will be scaled to fit with- 
in the specified frame size according to the 
following relationships. 

X-scale factor = (XRG-XLG) /SCAL 

Starting X-value = XLG 

Y-scale factor = (XRG-XLG) /SCAL 

Starting Y-value = XLG 

XRG and XLG are computed automatically by the IMAGE program. 

The coordinate axes used are the result of the initial geometric 
transformation described above. The relationships used for com- 
puting the above parameters internally are as follows: 

XLG = min (all X-coordinates) 

XRG = max (all X-coordinates) 

The above equation generates a "square frame" scaled to the 
longitudinal configuration geometry dimensions. The "square 
frame" produces identical scale factors in the X and Y directions. 
The geometry is guaranteed to fit in the specified frame size. 
However, the data will not necessarily fit in the frame itself. 

The position of the picture depends upon the values of the view- 
ing angles. For example, assume the geometry reference is at 
the nose of the vehicle. The imaginary frame under automatic 
scaling conditions would encompass the length of the vehicle. 
Further, the data would all be scaled to the value of the input 
variable, SCAL. For the above conditions, the position of the 
nose is controlled by the input variables DXG and DYG. 
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COORDINATE SYSTEM DEFINITIONS 

Xcc, Ycc Reference coordinates of the plotting device 

before the current picture is drawn. The 
initial values are usually specified by the 
user on the plot request card. Before the 
picture is drawn they are moved to the Xg, 

Yg system. 

Xg, Yg Reference axis of the geometric data after 

the initial transformation specified by the 
input parameters DELX, DELY, DELZ , XSC , YSC 
and ZSC. 


FIGURE 9 ILLUSTRATION OF THE FRAMING TECHNIQUE EMPLOYED IN 
THE IMAGE PROGRAM. 
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Each picture is positioned within the limits of the plotting 
device using the variables DXG and DYG measured with respect 
to the previous position. 

DXG = 0 . Repositioning of the reference Y-axis in 

inches for the plotting device before plotting 
the geometry data for the current picture. 

DYG = 0 Repositioning of the reference X-axis in 

inches for the plotting device before plotting 
the geometry data for the current picture . 

Once the reference axes for the current picture is described, 
the imaginary frame is defined by the variables YLG, XRG, YBG, 
YTG. Although computed automatically as described above, the 
values of XLG, YRG, YBG and YTG may be specified by the user 
in scaled geometry coordinates with respect to the translated 
vehicle geometry on TAPES as follows: 

XLG = computed Value of the left side of the imaginary 
frame (geometry scale) . 

XRG = computed Value of the right side of the imaginary 
frame (geometry scale) . 

YBG » computed Value of the bottom of the imaginary frame 
(geometry scale) . 

YTG = computed Value of the top of the imaginary frame 
(geometry scale) . 

When specifying the above values, the user should remember the 
"square frame" is essential to the creation of an undistorted 
picture. Rectangular frames will produce pictures which appear 
for-shortened in the direction of the short side of the 
rectangle. For example; 

XLG = -10. , XRG = 20. , 

YBG = -5. , YTG = 25. , 

would, produce a square frame 30 inches on a side. In the above 
example, the data for the current picture would be scaled as 
follows : 

Y - scale factor = (20 .- (-10 . ) ) /lO . = 3 units/inch 

Y - scale factor = (25.- (5. ) )/10. = 3 units/inch 
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starting X-value = -10 . 

Starting Y-value = -5. 

The primary function of the IMAGE program is the generation of 
pictorial data. Equally important to a geometric description 
is a description of the vehicle characteristics. The IMAGE 
program provides the user the option of displaying text informa- 
tion along with the pictorial data. The logical variable con- 
trolling this option is TEXT. 

TEXT = .FALSE. Logical variable, if .TRUE., Text informa- 
tion will be read from cards following the 
$TYP3435 data set in free field format. 

In the text option the input variables DXG and DYG position the 
reference coordinates at the beginning of the first line of 
text. This reference point will remain until altered by a new 
$TYP3435 data set. Any number of cards may be read as text. 

Each card represents a line of text. Lines of text may be 
skipped by placing a zero (0) in column 1 of the text card. 

This provides a convenient means of spacing the text information. 
The text input is terminated by placing the character (2) in 
column 1 of the last text card. The last card will not be 
printed. 

The maximum number of cards (lines) of text is 53 cards. The 
character height of the text is controlled by the input var- 
iable HTEXT. 

HTEXT = 0.14 Height of the characters in the text 
material . 


Post Processing Instructions 

The IMAGE program is designed to generate a file of plot commands 
for a 12-inch vertical height continuous roll paper such as 
CALCOMP. The file must be on a physical tape for plotting (see 
figure 3) . The user establishes a reference point on the roll 
at the time the plot request is submitted to be plotted. Gen- 
erally a one inch offset from the X-axis (Y=l) is adequate. 

The Y-offset is not applicable since a continuous roll of paper 
is generally used. A hard coded program header is printed 9.5 
inches above the user established reference. All pictures should 
be scaled and positioned within the 9.5 inch limitation on 
vertical height, otherwise, the picture may overlay the program 
header. The horizontal placement of pictures is essentially 
unlimited. 
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Externally Generated Geometry 

The IMAGE program provides the flexibility of reading the 
element data for plotting from the normal input device (INPUT) 
or from an alternate unit. The later capability is provided 
for the purpose of plotting geometry generated by another com- 
puter program. 

The information may be stored temporarily or permanently on a 
system file and attached to the IMAGE execution job step by a 
method called file substitution. 

File substitution is a GDC 6600 system capability which provides 
a correspondence between internal (logical units) files and 
external (system) files. The mechanism by which this correspond- 
ence is implemented is the "program card." The program card 
for the IMAGE program is : 

PROGRAM IMAGE (INPUT, OUTPUT, TAPES = INPUT, TAPE6 = 

OUTPUT, TAPES, TAPE 8) 

In the above illustration the file, TAPES is the internal 
logical unit which may contain the geometric data to be plotted. 
The correspondence between the internal logical units, TAPES 
and the external file is established by the substitution of 
the TAPES parameter at execution time. For example, assume 
the data to be plotted was stored on an external file called 
DATA. The execution card for the IMAGE program would be: 

EXECUTE (IMAGE, ,,,,,, ,DATA) 

During the above execution, the IMAGE program would read from 
the file called DATA each time the logical unit, TAPE 8, was 
read. The program card parameters are positional . Therefore, 
the six commas are essential for the proper use of file sub- 
stitution, one for each of the other files on the program card. 


Use of IMAGE within ODIN 

The Optimal Design Integration (ODIN) system is a library of 
independent computer programs representing the analytical cap- 
abilities in a wide variety of technological disciplines. The 
IMAGE computer program is but a single member of the ODIN library. 
The sequence of execution of the individual programs is con- 
trolled by the executive program, ODINEX (reference 2) which 
also maintains a name-oriented data base of design information. 
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Each piece of information is stored by name. The data base 
forms a communication link among the programs in the library. 
When used within the ODIN system, IMAGE receives data from 
the data base before execution. Generally, the ODIN library 
programs provided information to be stored in the data base. 
However, IMAGE does not currently generate information of this 
category. 

The actual transfer of information from the data base to IMAGE 
is performed by ODINEX through pre-processing of the data so 
the program is "unaware" that it is part of an analysis involv- 
ing many programs. There are no special input requirements for 
using IMAGE within the ODIN system. A single control directive, 

'EXECUTE IMAGE' 

is required for the execution of the program. The delimiter 
(') is a 4-8 punch. The data which follows this directive is 
the normal input data described above. However, any data 
values may come from the data base by specifying the data base 
name on the input card (in lieu of the actual value) . 

$TYPE32 


SCAL = ' SCALE ' , 


$END 

In the above illustration, the name SCALE is a data base name 
which may represent a scale factor for the entire configuration. 
The executive program > ODINEX, replaces the name and the 
associated delimiter 'SCALE' with the data base value for SCALE. 
Upon execution of IMAGE, the input component is photometrically 
scaled by the current scale factor in the data base. A pro- 
cedure is also available for transferring data base arrays (see 
reference 2) . 

The illustration above applies to namelist input, but the pro- 
cedure for extracting data base information is equally 
applicable to formatted input. The field width is specified 
to the position of the delimiters (') as illustrated below: 
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DEN is assumed to be the name of a data base variable. The 
value of DEN will be placed (by ODlNEX> on the input card in 
the most significant (E and F) format left justified in the 
specified field. If the data base variables were an integer, 
the value would be right justified in the field (I format). 


PROGRAM OUTPUT 


The output from the IMAGE program is both printed output and 
plotted output. Geometric characteristics of the quadrilateral 
elements are available at the user's option. This data is 
presented through the normal output channels. The plotted 
data is in the form of a plot tape. The actual plots are 
obtained by separate submission of a plot request card. Upon 
submission of the request, the plot tape is processed on the 
CALCOMP plotting hardware. Pictures of the vehicle at pre- 
selected viewing angles, as illustrated in figure 10, may be 
generated on the CALCOMP device . Input errors can then be 
corrected before the data is submitted to another ODIN program 
for the aerodynamic or other technology calculations. 
Alternately, annotated report quality pictures (see figure 1) 
may be generated by combining the picture and text options in 
the program. 

The IMAGE program can also be used to obtain a detailed print- 
out of the properties of each quadrilateral element of the 
vehicle as illustrated in figure 11. The normal output may 
also include the accumulated surface area as illustrated in 
figure 12 and the number of elements for each section of the 
vehicle. If no print options are specified, only the picture 
number is printed. 


ODIN Output 

The current IMAGE program generates no output for the ODIN 
data base. It has been used primarily for displaying geometric 
data to be used by other technology programs. 
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FIGURE lOB PICTURE OUTPUT FROM IMAGE. (CONTINUED) 
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FIGURE 11 ILLUSTRATION OF QUADRILATERAL ELEMENT PROPERTIES PRINTOUT FOR IMAGE. 
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FIGURE 12 ILLUSTRATION OF ACCUMULATED SURFACE AREA PRINTOUT FOR IMAGE, 


SAMPLE CASES 


Usually the most difficult aspect of using any computer pro- 
gram for the first time is mental inertia involved. During 
the initial learning period, the user gains the necessary 
confidence required to obtain useful results from the program. 
The learning period varies with the size and complexity of 
the particular computer program and the individual involved. 
The IMAGE program is small and therefore, the learning period 
should be short. 

This section presents four example problems designed to help 
the first-time user in overcoming the initial mental inertia. 
Once familiar with the program input, the user will find that 
providing data to the IMAGE program is not unlike providing 
data for a configuration layout. 


Three-View, Nose Right 

The first example is a three-view of a suborbital maneuvering 
vehicle. The input data for this example is illustrated in 
figure 12. Also illustrated in figure 13 is an inset drawing 
resulting from the input data shown. The drawing is reduced 
several fold from the original 27.9 x 43.2 cm (11 x 17 in) drawing. 

Each line of information in figure 13 represents a "card" of 
input data except as noted in the following discussion. The 
first card is the title card. This hollerith information 
appears just below the program heading as illustrated in the 
inset drawing. The information presented can be up to 59 
characters of the TITLE card. The next card is the entire 
$TYPE32 data set. Note that no data is entered in the data 
set for this example. Therefore, the program default values 
are used. By omitting the data entries in this section, the 
following input data is implied; 


PRINTS = 0, 

(no printing of quadrilateral 

characters) 

lORIEN = 0, 

(normal cross section 

input) 


ISTAT3 = 1, 

(one vehicle section) 



ITAPE = 0, 

(geometric input data 

in UNIT 

5) 

IREW8 = 0, 

(rewind UNIT 8 before 

reading 

- no effect 

XSC = 1.0, 

(X - scale factor) 
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FIGURE 15 ILLUSTRATION OF INPUT DATA FOR THREE-VIEW DRAWING - NOSE RIGHT 




YSC 


=1.0, (Y 
ZSC = 1.0, (Z 
DELX = 0 . , (X 
DELY = 0 . , (Y 
DELZ = 0. , (Z 


- scale factor) 

- scale factor) 

- translation) 

- translation) 

- translation) 


Even though no data is placed in the $TYPE32 data set, the 
empty data set shown must be present. The minimum information 
that must appear in the dummy data set includes the name, the 
opening $ and closing $ as follows: 


$TYPE32 $ 


At least one space must appear after the NAMELIST name. 

The next set of data is the geometric input data indicated by; 


GEOMETRY DATA 


The actual data is omitted from the figure for clarity of pre- 
sentation. Appendix A is a listing of the actual data used in 
this and other examples. 

The next series of data sets are the $TYP3435 data sets. These 
data control the position and orientation of a sequence of 
pictures representing the data described above. The text option 
is also illustrated. The correspondence between the $TYP3435 
data sets and the pictures (text) is indicated by the circled 
numbers. The crosshairs on the individual drawings refer to 
the location of the reference system at the start of the 
individual picture (text) . These locations are controlled by 
DXG and DYG. 


The first $TYP3435 data set defines a plan view orientation 
with the reference axes located 15.2 cm (6 in) to the right 
and 15.2 cm (6 in) above the initial ..reference system of the 
plotting device. The data is to be scaled from the longest 
X-dimension to fit within a 25.4 cm (10 in) imaginary frame. 
All other data in this data set will be the default values 
as follows: 

THETA = 0 

ICS = 0 


(pitch angle in degrees) 
(connect all four points) 
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IREFL = 0 
ISHAD = 0 
lAREA = 0 
IQUAD = 0 
LAST = 0 
XLG = computed 
XRG = computed 
YBG = computed 
YTG = computed 
TEXT = .FALSE. 
XMOVE =17.0 

YMOVE =0.0 
HTEXT =0.14 


(draw reflected plane) 

(omit drawing rear facing elements) 

(no section areas will be printed) 

(actual corner points will be drawn) 

(this is not the last picture) 

(left side of frame) 

(right side of frame) 

(bottom of frame) 

(top of frame) 

(this is a picture option, not text) 

(move plot device 17 inches in the X- 
direction after completion of the present 
picture series) 

(plot device does not move in the Y- 
direction after the current picture series) 

(character height for textual information) 


The second $TYP3435 data set defines a profile view of the 
vehicle. This view is to be located 12.7 cm (5 in) below the plan 
view (described above) . Note that all data remains unchanged 
between data set definitions. Only those variables requiring 
change from the previous set need be reset by the user. 

The third $TYP3435 data set defines a front view of the vehicle. 
This view is to be located 7.62 cm (3 in) to the right of the 
profile view. 

The fourth $TYP3435 data set defines a text option, a series 
of text cards which will be read from cards in 80 coliamn format. 
Only five input variables have meaning to the text option. 

TEXT (activates text option) 

DXG (X-movement of the plot device reference to the 

start of the first line of text) 
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DYG (Y-movement of the plot device reference to the 

start of the first line of text) 

LAST (indicates the last $TYP3435 data set in the 
current series) 

HTEXT (character height for textual information) 

The actual text immediately follows the $TYP3435 data set which 
activates the text option. The first column of each text card 
is reserved for print control as follows: 

0 - skip a line 

2 - terminate the text option 

In the illustration, the start of the first line of text is 
positioned 17.8 cm (7 in) aoove and 8.1 cm (2 in) to tne left 
of the previous (front view) reference and the characters are 
.36 cm (.14in). high. 

The parameter LAST is set to 1 in the text option indicating 
that current set to be the last $TYP3435 data set in the series. 
The program flow logic will return for a new title card. Then 
since the next (title) card has the characters "99" in columns 
71 and 72, the IMAGE program execution is terminated. 


Three-View, Nose Left 

The second example is a three view of a shuttle orbiter-type 
vehicle. The input data for this example is illustrated in 
figure 13. The drawing resulting from the data is inset into 
the figure. Except for the geometric data employed, the pri- 
mary difference between this example and the previous one is 
the orientation of the views. This example uses a nose left 
orientation and a rear view. The previous example used a nose 
right orientation and a front view. Comparison of figures 12 
and 13 will illustrate the differences between the data 
definition for the two types of vehicle orientation. 

The geometric data for this example was generated by a special 
computer program called PANEL (reference 3). Therefore, the 
$TYPE32 data set specifies the data source to be external by the 
parameter. 

ITAPE = 1 

This parameter specifies the data is to be read from the internal 
file, UNIT 8. The correspondence between UNIT 8 and the external 
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file which contains the actual data is established through the 
sixth parameter on the program card. The external source of 
data is accessed by the IMAGE execution control card through 
file substitution of the parameter with the external file 
name as follows : 

OIMAGE, , , , , ,GEOM. 

In the above illustration the file, GEOM, contains the geometric 
data for the IMAGE program. Of course, the geometric data must 
have been placed on the GEOM file by the PANEL program in a 
similar manner as that illustrated above. 

OPANEL, , , , , ,GEOM. 

In the case of the PANEL program, the geometry unit, GEOM, 
happens to be the fifth file parameter and must have been sub- 
stituted accordingly at the time PANEL was executed. 

The $TYPE32 data set also specifies a translation of the original 
data in GEOM (TAPES) 3048 cm (1200 in) forward and 1778 cm (700 
inches) down. The reason for the translation is to reposition 
the reference axis system for plotting purposes. The original 
data for the orbiter generated in the PANEL program was 
referenced to the coordinate system of the boost vehicle upon 
which the orbiter was mounted. The translated reference system 
is coincidental with the nose of the orbiter as indicated by 
the inset to figure 14. 

The picture sequence and text option which are defined by the 
$TYP3435 data sets are similar, except for orientation, to the 
sequence described in the previous example. Therefore, the 
reader is referred to that section for discussion of the input 
data. 


Oblique Views 

The third example shown in figure 15 illustrates a series of 
oblique views of the suborbital maneuvering vehicle of example 
1. The geometric data is read from the normal input device 
and not translated as was the data of example 2. The actual 
data is shown in Appendix A. 

Each $TYP3435 data set refers to a numbered view. The 
correspondence of the individual data set with the view is 
indicated in the figure. The significant difference between 
these data sets and the data sets of the other examples are: 
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ILLUSTRATION OF INPUT DATA FOR A 
SEQUENCE OF OBLIQUE VIEWS. 




1. Differences in viewing angles. 

2. Different locations for the views. 

3. The absence of text in this example. 

4. The frame dimensions XLG, XRG, YBG and YTG are read 
in rather than computed by the program. 

Notice the frame dimensions are established so that the frame 
size in the X-direction is the same as the frame size in the 
Y-direction. 

XRG - XLG = YTG - YBG 

The above criteria establishes a "square frame" which is 
essential to an undistorted picture. 


ODIN Input Example 

This example is an illustration of a three-view drawing of a 
shuttle orbiter recently studied at NASA Langley Research 
Center using the ODIN procedure to simulate certain aspects 
of the design process. Figure 16 is an illustration of the 
data setup for IMAGE within the ODIN framework. This figure 
portrays a set of data describing a three view, nose left with 
vehicle characteristics printed on the picture. The resulting 
picture is inset in the figure. Both the picture definition 
data and the optional text information is augmented by data 
base variable names. The variable names are data base inter- 
faces and therefore delimited by (^) . For details of the 
interface language, the reader is referred to reference 2. 

The primary differences between this example and the previous 
one illustrated in figure 13 are: 

1. The present example is an illustration of the program 
being executed within the ODIN system. The illustrated 
data was extracted directly from an ODIN simulation. 

2. The data is generated by the program of reference 1 
which was also executed in the ODIN simulation. 

3. The frame size is input from the data base and is based 
upon the current fuselage length (XLFUS from the data 
base) . 
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--^EXECUTE imaged 

^ADD DM=XLFUS/SCAL# _ . 

- ORBITER design 25K P/L _1.._ INCH=_«M . *. .. INCHES 

-STYPE32 ITAPE=1» ISTAT3=3« % 

STYP3A35 


i-PSI = -90.* 

-PHI = -90, » 

=.--- DXG=0., 


DYG=6., , _ 

===: XRG=/XLFUS ** 

-J XLG=0., 

^-1 YBG=0,, 

-J; YTG=)‘XLFUS Jt. 
-=SCAL = ^SCAL 
— S 

STYP3A35 
:^PH1 = 0,. 

_ DXG = 0.» 

_OYG = -5,, 

_S 

__$TYP3435 
-PSI = 180.* 

_ DXG=12,, 

OYG = 0,» 

-$ 

-STYP3A35 

--TEXT = .TRUE.* 

:^LAST" = 1, 

0XG=-2,5* 

__ DYG=7.5* 

_S 



_ VEHICLE characteristics 

- O '- . . 

- 1, MASS properties 

-»■ - - ^ 

landed WEIGHT *WLAND * 

entry weight - . i*wentRY^ 

- - DESIGN CG* SUBSONIC* 25K P/L *XCGTN * 

DESIGN CG hypersonic* 13K.P/I *XCGHYP*?._ 


- 4 . ... - - 

‘-,2, GEOMETRY 

.0 

AADO DM=XLFUS/I2.^ 

BODY length __ 

TOTAL WING AREA 
_ CHORD* THEO root 

CHORD* TIP 
ASPECT RATIO 

leading edge sweep angle 

TRAILING edge SWEEP ANGLE 
PADO DM=CRE*CTE*SPEXP/2.< 

ELEVON AREA* TOTAL 
EXPOSED wing location 


^DM * 
#ST0TAL^ 
WCROOT * 
*CJJP -* 
^AR . * 

aswtle * 

^SWTTE * 

Tom “ < 

#XOF , *_ 


*3,' AERODYNAMIC characteristics _ V . " 

0 . 

design trim lift coefficient ,^CLT(8)^ 

L DESIGN MINIMUM LANDING SPEED WVMO * 

MAX trim alpha hypersonic design CONO rfALTRiM# 


END OF PICTURE PHASE FIGURE 16 ILLUSTRATION OF INPUT 

DATA FOR imGE WHEN USED IN ODIN 
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4. The text option display actual vehicle characteristics 
from the data base for the current vehicle. The 
delimited names denote the data, base names for 

the indicated quantities . 

Each line of information in figure 16 represents a card of 
input data. The first card illustrated is the ODlisiKX control 
directive . 

T^EXECUTE IMAGE?^ 

The control directive is input information to the ODINEX 
executive system and accomplishes the following functions: 

1. Retrieves the IMAGE program from permanent storage at 
the beginning of the simulation. 

2. Executes the IMAGE program providing the necessary 
file substitution to pass geometry file to the program 
from the program of reference 1. 

3 . Provides for the return of control to the ODINEX 
executive system after completion of the execution of 
IMAGE. 

The second card is an ADD command which is part of the communica- 
tion language in the ODINEX executive system described in 
reference 2. The function of this card, 

7«ADD dm = XLFUS/SCAL?^ 

is to compute a scale of the drawing (DM) to be used on the 
picture title card. In the illustration the picture scale is 
1 to 132- units. The card illustrated above is "removed from 
the input stream by ODINEX and is therefore not read by the 
IMAGE program. 

The third card is the TITLE card. This card is generally the 
first card (when IMAGE is used as an independent program) . 
However, when used with the ODINEX executive system, the TITLE 
card may not be the first card due to interface requirements 
with- other programs and with the ODIN data base. However, after 
the input data is preprocessed by ODINEX, the modified input 
file will be identical in format to the previous examples. No 
delimited information will appear. 

The $TYPE32 data resets the parameter, 

ITAPE = 1 
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which specifies the geometric data will be obtained from an 
external source. The ODINEX executive system automatically 
maintains the proper correspondence between the IMAGE data 
file and the data files of the program of reference 1. The 
user need not be concerned with file substitution when using 
IMAGE with the ODIN system. 

Another parameter set in $TYPE32 is: 

ISTAT3=3 

This variable specifies that three sections of data on the 
geometry file will be plotted. Actually more data resided on 
the geometry file for the example simulation but the first 
three sections (WING, BODY and VERTICAL TAIL) were the only 
ones of interest for plotting purposes. 

The next series of data sets ($TYP3435) control the position 
and orientation of a sequence of pictures and text. The first 
set in the series also defines the frame size and dimensions. 
Frame size is referenced to the fuselage length, XLFUS and 
the data is scaled to the parameter, SCAL. Both of these 
parameters come from the data base as indicated in figure 16. 

The last $TYP3435 data set activates the text option. The 
actual text cards follow the $TYP3435 data set. These cards 
are modified by data base information as indicated by the 
delimited variable names. Each delimited name is replaced by 
the current value of the variable from data base. ODINEX 
performs this replacement function before the data is read by 
IMAGE. In this manner, the current vehicle characteristics 
are extracted from the data base and printed with the geometric 
representation of the vehicle. 
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APPENDIX A 


PART I 


LISTING OF TEXT DATA FOR THE IMAGE PROGRAM 


The following pages present a list of the element data for 
the suborbital maneuvering vehicle test data discussed in 
connection with the use of the IMAGE program. This sample 
data is also used in examples contained in Part II of the 
present report. 
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PART II - PROGRAM HIDDEN 


Program HIDDEN provides a true hidden line picture drawing 
capability for the ODIN system. This program differs from 
IMAGE both in computer time requirements and fidelity of the 
picture. HIDDEN requires significant amounts of computer 
time for execution but provides an accurate rendition of complex 
geometries. IMAGE, on the other hand, requires little computer 
time but only provides an accurate picture for convex objects. 

Thus, the combination of programs provides either accurate/ 
time consuming pictures or rough first cut/rapid running pictures. 
Choice' of programs in an ODIN simulation therefore rests with 
the analyst and his requirements. 

Program HIDDEN is an extended version of a program originally 
written at Purdue University. The present authors wish to 
acknowledge their debt to this source at this point. Use of 
existing codes is of course a major feature of the ODIN system. 

In fact, the ability to incorporate existing codes in the ODIN 
system is the underlying reason for the order of magnitude 
reduction in program costs when ODIN is compared to other major 
vehicle design synthesis programs which operate or seek to 
operate at the technology levels employed in current ODIN simu- 
lations . 

The description of Program HIDDEN which follows is an edited 
version of the Purdue University final report of Grant NGR-15- 
05-180 augmented by a description of the procedures used to 
generate pictures applicable to ODIN simulations. When used in 
the PCSYS picture drawing system HIDDEN operates with two addi- 
tional programs written by Aerophysics Research Corporation. 

These programs are CVRTHD which automatically converts IMAGE 
input to HIDDEN input and PLOTHD which converts HIDDEN output 
to a form acceptable to the Langley Research Center and Aerophysics 
Research Corporation plotting programs. Program CVRTHD and 
PLOTHD are described in Part III of the present report. 
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PROGRAM DESCRIPTION 


Coordinate System 


Program HIDDEN input consists of a list of points and a 
manner for specifying the connectivity of these points to 
form an object. All objects can be reasonably well repre- 
sented by lines, planes, and holes in planes. Thus, the most 
basic input scheme for a true hidden line algorithm requires 
a means of inputing points and then describing which points 
formed planes, holes, and lines. 

The specification of three dimensional points required the 
adoption of a reference coordinate system. The right handed 
coordinate system as shown in Figure 17 was adopted. It 
should be noted that this convention differs by a rotation from 
that of program IMAGE on page 3 . 



FIGURE 17. BASIC COORDINATE SYSTEM 
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The positive Y-axis can be thought of as being vertical and 
aimed upward and the X-Z plane to be the horizontal or 
ground plane. These statements are based on the viewing and 
display transformations present in the Hidden Line Removal 
Program. 


Surface Model Description 


When deciding the inputs to enter, a dominant goal was to 
keep the required data to a minimiam and to make the input as 
"natural" as possible. To make the input "natural" Purdue 
decided to use the words "POINTS", "LINES", "HOLES", "PLANES", 
and "END" to indicate the type of data following. These 
words are referred to as keywords. Each keyword sets the 
data type which is to follow and that data type remains in 
effect until another keyword is encountered. Thus, after 
the occurrence of the keyword "PLANES" several planes may 
be entered without each being preceded by the word "PLANES". 

If the planes, holes, and lines are to be described by 
references to points which have been previously read or are 
to be read, some scheme for the identification of the points 
must be included in the input routine. An obvious scheme is 
to identify each point with a number. The computer numbers 
the points sequentially, beginning at 1, as it encounters them. 
The user, although not required to enter these numbers, must 
be aware of them in order to enter the planes, holes, and lines. 
This solution avoids the problems arising out of the less 
ordered numbers a user might supply. It also saves some data 
preparation time. 

With each point assigned a sequence number all the user 
specifies a plane or hole by listing the sequence numbers 
assigned to the vertices forming the plane or hole. The sequence 
numbers must be entered in order proceeding around the surface 
element boundary. There is no restriction on the direction 
in which the boundary is traveled. Many algorithms require that 
the vertices be entered so that if vectors were formed from 
vertex 1 to 2, 2 to 3, etc., and two sequential vectors are 
crossed that an outward normal would be generated. This require- 
ment can be awkward for some three dimensional object input 
requirements. Program IMAGE avoids this problem by virture of 
other programs in the ODIN system which automatically prepare 
the surface geometries. Program PANEL, reference 6, is an 
example of one such program. In program HIDDEN the list of 
vertices forming planes or holes is completed with the second 
occurrence of the first vertex. Each new plane, hole, or line 
begines on a new line or card and, in the case of planes and 
holes, if there are more vertices than room permits the list 
is simply continued onto the next line or card. 
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If there exists other lines in the desired object besides 
those forming the plane and hole boundaries, their endpoint 
sequence numbers are given following the entry of a "LINE" 
keyward. The line feature permits the decoration or lettering 
of plane surfaces as well as the inclusion of such elements 
as ropes, wires, cables, etc. into three dimensional objects. 

Restrictions on data input are that holes must immediately 
follow the plane in which they occur and that after all data 
describing an object has been given an "END" keyword should 
be given. 


Data Formats 


All existing implementations of this basic input scheme have 
been done in FORTRAN- The format for the various pieces of 
data are as follows: 


KEYWORDS 

PLANES 

HOLES 

LINES 

POINTS 


A3,77X 
8X,1814 
8X,1814 
8X,214 
10X,3F10. 3 


An example is given in Figure 18 showing a simple object and 
the required data input. 


Program Operation 


As can be seen by the flow diagram of Figure 19, the program 
begins by asking the user to assign three files to logical 
unit numbers. Two of the three are for input and the remaining 
file is for output. One input file contains the object descrip- 
tion using either the basic input or the Composite Object Input 
Language. The other input file contains the viewing vector 
for each view. The viewing vector consists of two points given 
relative to the object coordinate system which denote where 
the observer is standing and where he is looking. The only 
restriction on this vector is that the observer must be outside 
the object because no clipping is provided. The third file, 
the output file, is that file which is written with the lines 
to be drawn. This output file can then be used as input to 
programs producing displays on suitable graphics display devices. 
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FIGURE 18. EXAMPLE OF BASIC INPUT 
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Subroutine INPUT2 


Subroutine INPUT2 contains the code for processing the object 
description as expressed using the basic input scheme. As 
the data is being read the coordinate arrays X, Y, and Z are 
filled as well as the line endpoint arrays LA and LB. The 
plane and hole vertices are stored on a temporary file. This 
is done to save core and is feasible because the vertices are 
only referenced a few times whereas the coordinates and line 
endpoints are continually being references during program 
operation. When the routine encounters an "END" keyword or 
the end of the object input file a return to the main program 
is executed. 


Subroutines INPUTl, COIL, and CNVRTl 


INPUTl is the subroutine which processes objects which have 
been described using the Composite Object Input Language. 

This mode of input has not been used in current ODIN appli- 
cations. The generalized data output arrays are filled in 
preparation for the CNVRTl subroutine. Subroutine COIL and 
the INPUTl were designed after the Hidden Line Removal Program 
to provide general three dimensional input and it was felt 
that the arrays that INPUTl produces should be more generalized 
than those required for the Hidden Line Removal Program. Thus, 
subroutine CNVRTl is required to convert the data formed by 
subroutine INPUTl into the required data. The required data 
being the file containing the vertices, the coordinate arrays, 
and the line endpoint arrays. 


Subroutine INTRSC - Plane Intersections 


Having processed all the object data the program generates 
all lines of the intersection formed by the intersections of 
planes. The intersection problem is handled in subroutine 
INTRSC. The general approach to the problem is to compute 
all lines of intersection plane 1 forms with plane 2 and sub- 
sequent planes, then to compute all lines of intersection plane 
2 forms with plane 3 and subsequent planes, and so forth. The 
addition of lines of intersection to the data simply involves 
the addition of new points to the coordinate arrays and new 
sequence numbers to the line endpoint arrays. The points are 
determined simply by computing all points where the boundaries 
of each plane pierce the other. If a boundary line does 
pierce the other plane then it is known that it is one point 
on a line of intersection. The analytics of a line piercing 
a plane are quite straight-forward and are discussed in 
Appendix B. Once all points along the intersection have been 
determined, the various pairs of points to be connected 
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1 



FOR THE SAKE OF THIS DIAGRAM CONSIDER THAT 
THE BOUNDARIES OF TWO PLANES. PLANE A AND 
PLANE B. HAVE BEEN FOUND TO PIERCE EACH 
OTHER NINT NUMBER OF TIMES ALONG THE LINE 
OF INTERSECTION. 


FIGURE 20. INTERSECTION LINE CONNECTIVITY LOGIC 


71 






becomes the problem. If a,ll the piercing points are placed 
in a list ordered by their position along the line of inter- 
section then the procedure shown in Figure 20 will generate 
the appropriate pairs. This logic is adequate enough to 
figure the connectivity of the lines of intersection which 
arise from the most complicated intersecting planes. 


Subroutine REDUND 


The lines which form the boundaries of the planes and holes 
must also be added to the line tables. This function is per- 
formed by subroutine REDUND. The subroutine received its 
name because of a constraint that is enforced. That constraint 
is to avoid adding redundant lines. When considering polyhedra 
it is obvious that most planes share boundaries, therefore, 
if the boundaries were added indiscriminately the number of 
lines to be processed and output would be almost doubled what 
needs to be. 

The coordinate and line endpoint arrays are then complete. 

These arrays will be altered then by transformations and the 
visibility process so if multiple view generation is to be 
permitted it is advantageous to save the data at this point. 
Thus, the coordiante arrays and line endpoint arrays are copied 
out to a scratch file. 


Viewing Vector and Viewing Plane 


The viewing vector is read next. This vector is used then to 
determine the coordinate transformation and perspective pro- 
jection for the view. The general purpose of the coordinate 
transformation is to translate the origin of the coordinate 
system to the point being observed and to rotate the coordinate 
system so that the positive Z axis contains the observer's 
position. The analytical method for performing this transform- 
ation is given in Appendix C. Once the coordinates have been 
transformed the perspective projection can be made. This 
essentially consists of projecting the object points onto a 
picture plane which is normal to the line of sight vector and 
passes through the point being observed. The resulting pro- 
jection is thus two dimensional and lies on the plane Z = 0. 

Two new arrays, XP and YP, are filled with the perspective 
representation of the coordinates. The actual analytical 
derivation and equations for the perspective projection are 
given in Appendix D. With the transformations out of the 
way the visibility problem is the next one to be solved. 
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Subroutine VISIBL 


The determination of which lines are visible or partially 
visible is perfoinned by subroutine VISIBL. The general 
approach taken in this subroutine is to consider each plane 
and determine which lines the plane hides or partially hides. 
The lines hidden or partially hidden are removed from the 
line endpoint arrays and if the line was only partially hidden 
the remaining visible segment or segments are added to the 
coordinate and endpoint arrays. After all the planes have 
been processed only those lines remain in the line endpoint 
table which should be displayed. Thus, the whole visibility 
problem reduces to the determination of methodology to resolve 
where a plane hides a line. The general outline for this 
methodology as implemented in subroutine VISIBL is shown 
schematically in Figures 21 and 22. 

As can be seen in Figure 21, the VISIBL subroutine begins 
with the code for referencing each plane. After a plane has 
been referenced the limiting coordinates are determined and 
an approximation for 0 is formed based on a percentage of the 
size of the plane. This approximation to 0 is required for 
subsequent calculations where round-off is a problem. The 
constants for the plane equation, AX + BY + CZ + D = 0, are 
then computed. Based on these constants the planes which 
appear as edge views may be trivially rejected because a 
plane which appears as an edge will not hide any lines . Next 
those planes which are part of convex objects and have had 
their vertices entered so as to generate an outer normal are 
considered. If the outer normal has a negative Z component 
the C in plane equation is negative and the whole plane need 
not be considered. This is so because the negative C means 
the plane is on the back of the convex polyhedron and anything 
the plane would hide will be hidden by the planes on the front 
of the polyhedron. There is no way to indicate a convex as 
opposed to a concave solid using the basic input scheme so 
this feature is only accessible through inputs such as COIL 
where the computer handles the computation of the direction 
of vertex input. 

Next the line endpoint arrays are compressed or have those 
lines eliminated which were indicated to be removed when the 
last plane was processed. The current number of lines is 
then stored because as the current plane is processed the 
total niimber of lines grwos with the generation of the line 
segments not hidden. There is then no need to consider these 
new lines against the plane when it is known they are not 
hidden by the current plane. Thus, by storing the number 
of lines in the table when the consideration of the plane 
begins it is knwon how many lines to consider even though the 
total number of lines increases. With the number of lines 
at the start stored, each line may then be considered. 
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FIGURE 21. SUBROUTINE VISIBL FLOWCHART 
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FIGURE 22. SUBROUTINE VISIBLE FLOW CHART (PART 2) 


75 















There are several trivial cases which can invmediately and 
easily be resolved. First, if the line is a boundary of the 
plane it cannot be hidden by the plane so the next line may 
be considered. If the line is a point view it need not be 
considered further and may be removed from the line table. 

Using the maximum and minimum plane coordinates determined 
earlier it is also possible to avoid further consideration 
of lines which do not even lie in the same reion as the plane. 

If a line is still to be considered after all of these constraints, 
more involved computations are required. 

The initial computations are involved in determining if either 
endpoint is behind the plane at least in depth. This first 
level of consideration ignores whether the line lies outside 
or inside the boundaries of the plane as projected on the 
picture plane and considers the problem as if a line and 
infinite plane were the case. The depth consideration is based 
on the distances from the observer to the line endpoints and 
to the projections of the line endpoints and to the projections 
of the line endpoints onto the plane. Thus, the first calcula- 
tions compute the distances from the observer to line endpoints 
1 and 2. Then it is attempted to project the points back to 
the plane. If the line through the observation point and line 
endpoint could not be projected onto the plane and is on the 
visible side of the plane. If the projection line is not 
parallel to the plane then the piercing point of the projection 
line and plane is computed along with the distance from the 
observer to this piercing point. If both line endpoints are 
in front of or on the plane, the plane cannot hide the line and 
the next line may be considered. If the line intersects the 
plane or lies entirely behind the plane, additional processing 
is necessary. 

If the line intersects the plane, point 1 will be behind the 
plane and point 2 will be in front of the plane. The piercing 
point is computed and then it is needed to know whether the 
piercing point occurs inside or outside the plane's boundaries. 

A piercing point which occurs within a hole or notch in the 
plane must obviously be considered outside the plane. 


Subroutine INSIDE 


Subroutine INSIDE decides whether the first point of the line 
supplied to it is inside or outside the plane's boundaries 
as well as returning all apparent intersection points between 
the line and plane boundaries as projected in the picture plane. 
INSIDE also indicates whether the line is in front of or 
behind the plane at each intersection point. If the piercing 
point falls within the plane boundaries, the original line is 
deleted and a line from the peircing point to endpoint 2 is 
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added to the line endpoint arrays with the piercing point 
being added to the coordinate arrays- Point 2 of the 
original line is replaced by the piercing point coordinates and 
the modified line is then considered as a line lying totally 
behind the plane as is described in the next paragraph. If, 
however, the piercing point falls outside the boundaries of 
the plane the line is sent to INSIDE to determine if point 1 
is inside or outside the boundaries of the plane. 

Again as well as resolving the intersection problem all apparent 
intersections of the line and plane boundaries along with depth 
information are returned. When the list of apparent inter- 
section points is returned line endpoint 1 is the first in 
the list, line endpoint 2 is the last in the list and all others 
are appropriately ordered in between. If point 1 was inside 
the plane's boundaries, lines are added from the even apparent 
intersections to the odd apparent intersections until the 
apparent intersections are in front of the plane. If point 1 
was outside the plane's boundaries, lines are added from the 
odd apparent intersections to the even apparent intersections 
until the apparent intersections are in front of the plane. In 
both cases when an apparent intersection was finally found which 
was in front of the plane a line was added from the previous 
apparent intersection to endpoint 2. 

When the line lies entirely behind the plane the procedure is 
greatly simplified over that required when the line intersects 
teh plane. The line which lies entirely behind the plane is 
first sent to INSIDE to find out if point 1 is inside or out- 
side the plane's boundaries and to determine the apparent 
intersections. If point 1 is inside the plane, lines are added 
from even apparent intersection points to odd apparent inter- 
section points until all the apparent intersections have been 
considered. However, if point 1 is outside the plane, lines 
are added from odd apparent intersection points to even apparent 
intersection points until all the apparent intersections have 
been considered. 

Whenever a newly generated line segment is added the original 
line from which the line segment was generated is deleted. 

The logic to decide whether the first point of a line is inside 
or outside the plane is much simpler than might originally be 
anticipated. The line is extended from point 1 through point 2 
until the new endpoint is beyond the plane. The maxima and 
minima computed earlier for the plane are employed for this. 

Then all apparent intersections are computed for the extended 
line and the boundary lines for the plane as they appear on 
the picture plane. Any hole boundary lines are also considered 
at this point because they too are boundaries of the plane. 

The apparent intersectionsare projected back to the plane and 
line and the distances are computed from the observer to the 
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points on the plane and line respectively. Then if the number 
of apparent intersections is odd the first point is inside 
the plane and if it is even the first point is outside the 
plane. 


Subroutines SIZE and DRAW 


After VISIBL has produced line endpoint arrays with only the 
lines to be displayed, subroutine SIZE can use this data 
to compute the maxima and minima picture plane coordinates 
of the points to be displayed. The maxima and minima can 
then be used along with the plot size to compute a scale 
factor. The minima and scale factor are then passed through 
common and argument list to subroutine DRAW. Subroutine 
DRAW writes out the X minimtim, Y minimum, and scale factor 
in a 2H 1,3F7.2 field and then outputs the lines. The lines 
are output by using the sequence numbers from the line end- 
point arrays to reference the pairs of coordinates from the 
perspective coordinate arrays. The pairs of coordinates are 
then simply written out in a 2X,4F.7.2 format. 

Having written out the information for the view, the original 
coordinate and line data is retrieved from the scratch file 
and the next view is considered. Upon completion of all 
views the program execution is terminated. 

The output file which was created can then be used as input 
by programs which drive the appropriate graphics device, for 
example, program PLOTTR, reference 7. These programs read 
teh scaling and minima information and use it in the plotting 
of the lines. 
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APPENDIX B - PART II 


COMPUTATION OF THE P IERCING POI N T OF A LI N E THROUGH A PLANE 


The computation of the piercing point of a line through a 
plane is nothing more than the simultaneous solution of the 
equations for the line and a plane. The equation of a plane 
is given by 


Ax + Ry + Cz + D = 0 

and that of a line by 

x-Xi ^ yrli_ = 2-Zi 
X 2 X 1 Ya'Yi Z 2 ~ Z 1 


( 1 ) 

( 2 ) 


where X^^, and X 2 , ^2 

Solving equation 2 for y and z in terms 


line endpoints, 
of X yields 


(x-Xi) 

X2-X1 

( 3 ) 

(x-Xi) 

+ Zi 

X2-Z1 

( 4 ) 


Substituting y and z as given in equations 3 and 4 into 
equation 1 yields 


Ax + B[(x-X.) + Y,] + C[(x-Xi) 

+ D = 0 
or 


x[A + 




BX, ( 


Y2-Y 

X2-X 


1 

1 


) 


BYi + 


CX, ( 



CZi 


D . 
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Solving for x 


X 


A + B(P^-p-) + C( 

A2 •" A 1 


Z 2 -Z 

X 2 -X 

T-PT 

X2-X 


1 

i 

1 


) 

) 



Then letting 


X2-X1 = V, 


Y 2 -Y 

X 2 -X 


i 

1 



Z 2 -Z 1 

X2-X1 



(5) 

( 6 ) 

(7) 


the equation for x becomes 

X = B(XiU-Yi ) + C(XiW-Zi ) - D . 

A~+ BU + CW 

Equations 3 and 4 may then be used to compute the y and z 
coordinate of the piercing point as 


y = (x-X 1 ) U + Yi 

and 


(9) 


z = (x-Xi ) W + Zi 


( 10 ) 


If V = 0, equations 6 and 7 yield «>. So the line equations 
are solved for x and z in terms of y and these equations are 
substituted back into equation 1. Equation 1 is solved for 
y and with 

V = Y 2 -Y 1 

B ” jLzT-Ai. = X2 - X 1 _ _ 

Y2-Y1 \r^ - 0 

W = ^ -_Zi _ Iz-Ii 

Y2-Y1 r~ 

yields 

y =- AXi + CTYiW- Zi] - D 
B + CW' 
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X = Xi or Xj 

z = (y-Yi ) W + Zj 


If both X 2 "X^ = 0 and ^ 2~^1 ° 


z = -AXi- BY,-D 
C 

and 

X = X 1 or X2 
y = Yi or Yz . 



APPENDIX C - PART II 


COORDINATE TRANSFORMATION 


In order to create the various views desired, the coordinates 
of an object must be transformed. The transformation consists 
of changing the coordinate system from its original orienta- 
tion to an orientation in which the origin is the point at 
which the observer is looking, (XL,YL,ZL), the positive Z 
axis is aimed toward the point where the observer is located, 
(X0,Y0,Z0). This transformation is thus a translation and 
two rotations as shown in Figure C-1. 



FIGURE C-1. COORDINATE TRANSFORMATION 
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The transformations to perform each of these tasks will be 
defined and then concatenated. 

The translation is the simplest transformation and can be 
expressed as 

X-XL = X' 

Y-YL = Y' 

Z-ZL = Z' 

in matrix form with homogeneous coordinates the transformation 
becomes 


1 0 0 -XL~ 


X 


X' 

0 1 0 -YL 


Y 

__ 

Y ’ 

0 0 1 -ZL 


Z 


Z' 

0 0 0 1 


1 


1 


or 

[Tj] [R] = [R*] 

Two rotations are required as shown in Figure C-1. The first 
is a rotation about Y' to bring Z' under the observation vector 
and the second is a rotation about X" to make Z" and the 
observation vector coicident. 

The problem is to determine the new coordinates (R') in the 
rotated coordinate system from the coordinates (R) in the old 
system. The general transformation, (Hall, 1966) , is given 
below for the rotation of the coordinate system about an 
arbitrary axis. 


(Uj^NersY + cosy) (Uj^UyVersy + U^si ny) vers y-UyS i ny ) 

(Uj^UyVerSY-U^siny) (U^yVersy + cosy) (UyU^versy+Uj^siny) 
(U^Uj^versy+UySiny) (UyU2.versy-Uj^si ny) versy+cosy ) 



[T] [R] 


[R"‘] 



where U is unit vector along the axis of rotation, y the angle 
of rotation, and vers = 1-cosy. 


If the transformation and vectors are extended to homogeneous 
coordinates : 


(Ujj^versy+cosy) ( Uj^UyVersy + U^s i ny) (Uj^U^versY-UySiny) 0 
(Uj^Uyversy-U^siny) ( Uy " ve rsy+cos y ) (UyU^versy+U^^siny) 0 
^ ^ ( Uy ve rsy- L) i ny ) { ve rsy+cos y ) 0 


•[T 


X 


’x'" 

Y 

/ 

Y' 


[R] = 


Z 


Z' 



1 _ 


[R] = 


Now if a rotation about Y is desired with an angle of rotation 
a. Then Uj^=U 2 = 0 and U^=l, y=a, the transformation (T) becomes: 


cosa 

0 

s i na 
0 


0 

1 

0 

0 


- s 1 n a 
0 

cosa 

0 


0 

0 

0 

1 


= [TJ 


If a rotation about X is desired with an angle of rotation 3. 
Then Uy=U 2 0, Y=3, and the transformation (T) becomes; 

To 0 0 

0 cos 3 sinp 0 

n = I^Ty] 

0 -sinp cos6 0 ^ 

0 0 0 1 

_J 
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Thus the total transformation can now be written as 
[Tj(] [Ty] [Ty] [R] = [R-] 

Expanding (T,p) (R) = 

X-XL 

Y-YL 

Z-ZL 

1 

Multiply that result by (T^) 

cosa(X-XL) - (Z-ZL) sina 
Y-YL 

(X-XL) sina + (Z-ZL) cosa 
1 _ 

Multiplying that result by (T ) 


(X-XDcosa-(Z-ZL)sina 


X'" 

(Y-YL) cose+( X-XL) si nab in6+ (Z-ZL) cos asinB 


Y *'* 

-(Y-YL)sinB+(X-XL)sin<.. cosp+(Z-ZL)cosucos8 


Z». 

_1 


J 
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APPENDIX D - PART II 


PERSPECTIVE TRANSFORMATION 


Once the Z axis is aligned with the line of sight vector the 
points must then be projected onto the YX plane to generate 
teh desired perspective view. This projection is shown in 
Figure D-1. 

The analytical derivation is quite straight-forward. Noting 
similar triangles in Figure D-1 the following equations may 
be written 


_YP Y 

DIST ■ DIST-Z 


and 


V p Y 

TOT " msT-Z 


Thus 


_ XDIST 
DIST-Z 


and 


_ YDIST 
DIST-Z 


86 
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FIGURE D-1. PERSPECTIVE TRANSFORMATION 
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PART III 


INTERFACING PROGRAMS IMAGE AND HIDDEN 


Aerophysics Research Corporation has constructed a program 
interface between the IMAGE and HIDDEN codes. This interface 
permits program HIDDEN to operate on program IMAGE data. 

Thus, when the ODIN system creates an input file for IMAGE 
the same file may be used to draw accurate hidden line 
pictures through HIDDEN. The two picture drawing codes 
are therefore available to ODIN system users with a unified 
input procedure. In view of the time consuming calculations 
required when producing true hidden line pictures through 
program HIDDEN a user has the option of "thinning out" the 
IMAGE data if desired. 
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IMAGE SUMMARY 


Image Input 

IMAGE input consists of corner points in the form 



where X is a coordinate normal to the paper plane. The hidden 
line algorithm uses Z as the normal to the paper plane, hence, 
CVRTHD assumes the data (on UNIT 5) has the fonr: 



The integer Ij^ determines corner point status as follows: 


I^ = 2, Start of new vehicle section 
= 1, Start of new row 
= 0 , Internal point 


An end-of-file test is used to define the end of data. It 
is assumed that the vehicle geometry is given only on one 
side of the plane of symetry. 


Program CVRTHD - Data Conversion Program 


Main program for converting IMAGE data to HIDDEN data. It 
uses two subroutines, PNTS and PLANE. The subroutine PNTS 
first converts the IMAGE corner points into HIDDEN format. 
Subroutine PLANE connects these points together to supply 
the "PLANES" data required by HIDDEN. A flow chart for data 
conversion is given in Figure 23. Data input to CVRTHD is: 

First Card (Format 12) 


IDIV = Averaging factor for reduction in number of points . 

End points are fixed but internal points are averaged 
by combining IDIV cards. Where the number of cards 
in a row leaves a remainder less than IDIV the 
program averages as many as possible . Some typical 
examples of this averaging process are shown in Figure 24. 
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Remaining Cards (Format 3(F9.3, IX), 1, 3(F9.3, IX), (II) 


The variables 



which successively define the panelled surface. 
Subroutine PNTS reproduces them in the form 


^1^ 

Xi, 

1— 1 

^2' 

X2, 

^2' 

Z 3 / 

X3, 

1 

^ 3 ' 

z , 

1 

X , 

Y 

n 

n 

n 


Subroutine PLANE then forms the triangular or quadrilateral 
panels using Z^, X^, i = 1, as illustrated in 

Figure 25. Planes are output as integers, for example. 



where these integers correspond to those on Figure 25b. Both 
points and planes data are output on file "CVTOUT." In 
Figure 25 (a) a typical panelling is illustrated with all 
points distinct. Figure 25(b) illustrates a case where one 
row of points are all equal. For example at a fuselage nose. 
A series of triangular panels are required in HIDDEN at the 
aircraft nose for zero length lines are not acceptable to 
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,-Unit 5 , 

I IMAGE INPUT I 


UNIT "CVTOUT" 
=6 in CVRTHD 
=lin HIDDNI 


Tape 4= INPUT 



Viewing 

Data 


1 


CVRTHD 
CONVERTS 
TO HIDDE 

IMAGE DATA 
DATA 


< 

FILE OF "POINTS 
AND "PLANES WHICH 
USE THE POINTS 



HIDDEN 

MAIN PROGRAM HIDDEN 
LINE ALGORITHM 

Tape 5 > 


1 PLOT FILE HIDOUT | 



PLOTHD 

PREPARES PLOT TAPE 
ON UNIT 2 


^ PNTS I 
PLANE I 


•Points 

Prepared 

•Planes 

Prepared 


S For each 
section o£ 
the vehicle 


•PURDUE 


■> 


PLOTTER 

PACKAGE 


FIGURE 23. SCHEMATIC OF HIDDEN LINE OPTION 
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FIGURE 24. TYPICAL 
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PLANES 


13 2 1 

14 3 1 

15 4 1 

16 5 1 

17 6 1 

18 7 1 

19 8 1 

2 3 11 10 2 

3 4 12 11 3 


I 


32 33 41 40 3; 


FIGURE 25(b). PLANE CONVENTION SEVERAL POINTS COINCIDENT 






HIDDEN. Figure 26 illustrates other simple panellings that 
may be encountered. Program CVRTHD automatically redefines 
quadrilateral panels as triangular panels wherever this is 
necessary. 


Program HIDDEN Input and Output 


Input 

The input to HIDDEN consists of points and planes on file CVTOUT 
as discussed previously. File CVTOUT ends with the keyword 
END starting in column 1. The format of this file is therefore 


POINTS 

Point triplets for 1st section 
PLANES 

Plane integers for 1st section 
POINTS 

Point triplets for 2nd section 
PLANES 

Plane integers for 2nd section 
I 
( 

( 

etc , 

POINTS 

Points triplets for last section 
PLANES 

Plane integers for last section 
END 


The second input file to HIDDEN is the last card of the input 
data deck. This card has the format 6F10.0. The six quantities 
on this card are 


XVi, YV^, ZV^, XV^, YV^, ZV^ 


where these two triplets define the viewing angle. The point 
XV^, AVj^ should be placed near the vehicle. The point 
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^2' ^2' ^^2 pla,ced at the viewers location. 

Data deck input for a hidden line run therefore has the 
form illustrated in Figure 27. The viewing vector card 
is internally transferred to Unit 4 in HIDDEN by file 
substitution. 


Output 

Program HIDDEN outputs the file HIDOUT on Tape 5. This file 
is then read by program PLOTHD in preparing the plot file. 
The contents of file HIDOUT are in the form: 

Record 1 - Scaling Data 

I J, XOFFSET, YOFFSET, SCa£e~| 


where , 

J = Record status flag 

= 1, First record-scaling data 
= 0 , Any other record 

XOFFSET = Internally computed offset for X 

YOFFSET = Internally computed offset for Y 

SCALE = Internal scale factor for adjusting 
picture size in HIDDEN (picture size 
can be modified in PLOTHD) 

Succeeding Records - Line Segments 



where , 

J = 0 

Xj^,Y^ are 1st endpoint of line segment 
^ 2''^ 2 endpoint of line segment 
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FIGURE 27. INPUT DATA DECK TO CONVERT IMAGE DATA TO HIDDEN LINE PICTURE 








Program PLOTHD - Picture Drawer 


Program PLOTHD reads the scale info 2 rmation and line segment 
corner points on file HIDOUT and prepares the X-Y plotter 
file using the Aerophysics Research Corporation Plotting 
Software Package (PLOTTR) , reference 7, or Langley Plotting 
Software. A flow chart of PLOTHD is given in Figure 28. 
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FIGURE 28. FLOW CHART FOR PLOTHD 
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APPLICATION TO THE FDL-6 CONFIGURATION 


AN ILLUSTRATIVE PICTURE SEQUENCE OF HIDDEN LINE PICTURES 


As an illustrative example of hidden line drawing ability the 
FDL-6 vehicle data used in Part I was supplied to PCSYS. 

The data was converted to HIDDEN format by CVRTHD. Program 
HIDDEN then prepared the hidden line data which was drawn 
on the Houston plotting device connected to a CDC 7600. The 
PLOTTR software was used to draw the resulting picture 
sequence for a series of viewing angles around the FDL-6 body. 
The fidelity of the resulting pictures is apparent in 
Figure 29. These pictures should be compared to those pre- 
sented in Part I, Figures 10(a) and 10(b), which were produced 
by the more rapid running program IMAGE. In particular, 
attention is drawn to the upper three-quarters rear view of 
Figure 10 (b) where the vehicle lower surface and interior 
is lost in the IMAGE picture. In the PCSYS pictures it is 
possible to see inside the vehicle from a reward direction 
since base panels were omitted. 
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FIGURE 29(i-l). PDL-6 CONFIGURATION 



105 







FIGURE 29 (m-p) . FDL-6 CONFIGURATION 


106 


/ 




FIGURE 29(q-t). FDL-6 CONFIGURATION 




APPLICATION TO AN ADVANCED TRANSPORTATION SYSTEM 


CONFIGURATION 


A final demonstration applies to IMSYS to an Advanced 
Transportation System configuration derived from the ODIN 
system by NASA personnel at Langley Research Center. The 
pictures obtained are presented in Figure 30. These pictures 
are notable for the ability of the codes to draw a complex 
picture illustrating four aileron positions and again to 
examine the vehicle interior from rearward positions. Once 
again these pictures were drawn from data originally supplied 
to the IMAGE code. Typical pictures provided by IMAGE from 
the same data are presented in Figure 31. 
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FIGURE 30(a). SOME TYPICAL ADVANCED TRi\NS?ORTATION VEHICLE VIEWS 



FIGURE 30(b). SOME TYPICAL ADVANCED TRANSPORTATION VEHICLE VIEWS 

FROM PROGRAM HIDDEN 
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FIGURE 31. CONFIGURATION FROM IMAGE 



CONCLUSION 


An improved picture drawing code (PC5YS) for the ODIN Optimal 
Design Integration System has been developed. The picture 
drawing code retains the input format of the original ODIN 
picture drawing program IMAGE which, in general, has been well 
accepted by users of the ODIN system. The improved picture 
drawing code has the ability to produce plots with no hidden 
line removal, with a rapid approximate hidden line removal algo- 
rithm based on inclination of the outward normal, or finally 
with an accurate but time consuming hidden line removal algorithm. 

This combination of picture drawing options enhances the graphics 
capability of the ODIN system while maintaining ease of user 
operation. The study has therefore contributed to the long 
range goal of improving and expanding the ODIN system technology 
modules. 
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